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ABSTRACT 

The intergalactic magnetic field (IGMF) may leave an imprint on the angular 
anisotropy of the extragalactic gamma-ray background through its effect on electro- 
magnetic cascades triggered by interactions between very high energy photons and the 
extragalactic background light. A strong IGMF will deflect secondary particles pro- 
duced in these cascades and will thus tend to isotropize lower energy cascade photons, 
thereby inducing a modulation in the anisotropy energy spectrum of the gamma-ray 
background. Here we present a simple, proof-of-concept calculation of the magnitude 
of this effect and demonstrate that current Fermi data already seem to prefer non- 
negligible IGMF values. The anisotropy energy spectrum of the Fermi gamma-ray 
background could thus be used as a probe of the IGMF strength. 

Key words: magnetic fields - galaxies: active - BL Lacertae objects: general - 

galaxies: starburst - diffuse radiation. 


1 INTRODUCTION 

Cosmic magnetic fields are expected to play a fundamental 
role in the physics of a large variety of astrophysic.al systems. 
The fate of high energy particles, such as ultra-high energy 
cosmic rays (UHECRs), as they propagate through the Uni- 
verse and the acceleration of charged particles in astrophys- 
ical objects hinge, in large part, on magnetic fields within 
their sources and those permeating the Universe. Large-scale 
magnetic fields such as those found in galaxies and clus- 
ters of galaxies have been a mystery for the past several 
decades. While they are generally thought to be the result 
of the amplification of weak seed fields, the origins of these 
seed fields, whether they be cosmological or astrophysical 
in nature, are largely unknown. A definitive measurement 
of the intergalactic magnetic field (IGMF) could provide a 
fundamental step in resolving the questions of the origins of 
cosmic magnetic fields and their impact on the evolution of 
the systems in which they reside, but sufficiently constrain- 
ing observations have thus far remained elusive. Previously, 
upper bounds of ~ 1CU 9 G have been found through Fara- 
day rotation limits of polarized radio emission from distant 
quasars and the study of the effect of magnetic fields on 
the anisotropy of the cosmic microwave background (for a 
brief review of observational and theoretical techniques to 
constrain the IGMF, see e.g., Neronov & Semikoz 2009, and 
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references therein). The recent availability of data from the 
Fermi Large Area Telescope (LAT) has renewed interest in 
the IGMF as observations of the gamma-ray sky and its par- 
ticipating source populations provide a unique perspective 
on the strength and structure of the IGMF. 

The gamma-ray sky consists of resolved point sources 
(such as normal and active galaxies, pulsars, etc.), transient 
gamma-ray sources (e.g., gamma-ray bursts), and the diffuse 
gamma-ray radiation comprised of emission from the Galaxy 
and the isotropic (presumably, extragalactic) gamma-ray 
background (EGB). The EGB is a window into the high- 
energy processes in the Universe, and its origins have been 
the subject of much debate. It is expected that emission 
arising from unresolved, extragalactic point sources, such 
as blazars and star-forming galaxies (SF), comprises a siz- 
able contribution to the EGB (see e.g., Stecker & Salamon 
1996; Fields et al. 2010; Abdo et al. 2010e; Stecker & Ven- 
ters 2011, and references therein). Additionally, many of 
these extragalactic point sources are also sources of very 
high energy (VHE) gamma rays, 1 which interact with the 
soft photons of the extragalactic background light (EBL), 
consisting of the cosmic microwave background (CMB) and 
the background of infrared, optical, and ultraviolet radiation 
from direct and dust-reprocessed starlight. The interactions 
of VHE photons with the EBL initiate electromagnetic (EM) 


1 In this paper, we take VHE to be ~ TeV. 
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cascades, giving rise to another contribution to the EGB (see 
e.g., Coppi & Aharonian 1997; Kneiske & Mannheim 2008; 
Inoue & Totani 2009; Venters 2010). However, the degree to 
which astrophysical sources contribute to the EGB depends 
on their distribution with respect to luminosity and redshift 
(the gamma-ray luminosity function, GLF), their intrinsic 
spectra at gamma-ray energies, and the nature of the EBL, 
all of which remain the subjects of intense debate (see e.g., 
Inoue & Totani 2009; Venters 2010; Fields et al. 2010; Abdo 
et al. 2010b, e; Makiya et al. 2011; Inoue 2011; Stecker & 
Venters 2011; Abazajian et al. 2011). 

A promising complementary technique for investigat- 
ing the contributions of astrophysical sources to the EGB 
involves studying the angular fluctuations in intensity of 
the EGB as a function of energy, the anisotropy en- 
ergy spectrum. As pointed out in Siegal-Gaskins & Pavli- 
dou (2009; see also Hensley et al. 2010; Siegal-Gaskins et al. 
2011), the contribution of a given class of sources to the 
overall fluctuation angular power at a given angular scale 
of the EGB at a given energy is weighted by the square of 
its fractional contribution to the intensity of the EGB at 
that energy. As such, if the relative contributions of the var- 
ious EGB contributors change with energy, so too will the 
fluctuation angular power as a function of energy. Expec- 
tations for the anisotropies of the commonly invoked astro- 
physical contributors to the EGB are that SF galaxies con- 
tribute very little anisotropy owing to the fact that they are 
so numerous and individually faint, while blazars, being far 
more rare than SF galaxies but much brighter, could have a 
large contribution to the EGB anisotropy (Ando et al. 2007; 
Ando & Pavlidou 2009; Siegal-Gaskins & Pavlidou 2009; 
Hensley et al. 2010; Ackermann et al. 2012a). It has also 
been suggested that annihilation in dark matter substruc- 
tures in the Galactic halo could also have a large contribu- 
tion to the EGB anisotropy, even if the contribution to the 
EGB intensity from dark matter annihilation is subdomi- 
nant with respect to other contributions (see e.g., Ando et al. 
2007; Siegal-Gaskins 2008; Taoso et al. 2009; Ando 2009; 
Siegal-Gaskins & Pavlidou 2009; Hensley et al. 2010). Early 
anisotropy results from the Eerrm-LAT Collaboration (Ack- 
ermann et al. 2012a) suggest that the energy dependence of 
the fluctuation angular power of the EGB may be consistent 
with arising from a source population with spectra similar 
to that of blazars (a ~ 2.4), even though blazars may com- 
prise < 30% of the EGB intensity (Abdo et al. 2010e; Cuoco 
et al. 2012). Thus, even contributions to the EGB that are 
subdominant with respect to other contributions can have a 
profound impact on the anisotropy of the EGB, especially if 
the relative prominence of its components changes with en- 
ergy. As such, the anisotropy energy spectrum could prove 
a powerful tool in detecting changes in the EGB anisotropy 
resulting from the magnetic deflection of EM cascades from 
blazars. 

In interacting with EBL photons, VHE gamma rays ini- 
tiate EM cascades of photons, electrons, and positrons. The 
interaction of a VHE gamma ray with an EBL photon pro- 
duces a pair of electrons and positrons, which will, in turn, 
Inverse Compton (IC) scatter EBL photons to high ener- 
gies. These up-scattered photons will, in turn, pair produce, 
and the process continues until the energies of the resulting 
photons are low enough that pair production is no longer 
efficient. In the presence of an IGMF, the charged parti- 


cles in the cascade are deflected away from the direction of 
propagation of the primary photon. As these charged par- 
ticles propagate through the EBL, they can up-scatter soft 
photons towards the observer, though they will appear to be 
coming from a different direction from the GeV gamma rays. 
Observationally, this process results in a “halo” of gamma 
rays around the source. The detection of a halo around a 
gamma-ray point source would not only provide strong evi- 
dence for an IGMF, but would also provide an indication of 
the strength of the IGMF. Several groups have already con- 
ducted searches for extended emission around gamma-ray 
sources (see e.g., Aharonian et al. 2001; Neronov & Vovk 
2010; Ando & Kusenko 2010; Aleksic et al. 2010; Neronov 
et al. 2011; Arlen et al. 2012), though as yet, no significant 
indication of gamma-ray halos distinguishable from instru- 
mental effects have been reported. However, the absence of 
significant halo emission could imply a strong IGMF that 
would extend the EM cascade halo such that the surface 
brightness of the extended source falls below the detection 
threshold (Neronov & Vovk 2010; Ando & Kusenko 2010). 
In addition to possibly creating a halo, the deflection of EM 
cascades by the IGMF could result in a reduction in ob- 
servable cascade radiation. Analyses employing simultane- 
ous observations of blazars by atmospheric Cherenkov tele- 
scopes and Fermi make use of this fact to derive limits on the 
IGMF (Dermer et al. 2011; Essey et al. 2011; Taylor et al. 
2011; Kachelriefi et al. 2012). These techniques are quite 
promising and, with more data from Fermi , will continue to 
shed light on the IGMF. In this paper, we consider a dif- 
ferent promising technique that makes use of the anisotropy 
energy spectrum of the EGB to constrain the IGMF. 

For any cosmological population, such as blazars, that 
emits gamma rays at very high energies, the effect of EM 
cascading results in a flux suppression at the highest en- 
ergies and enhancement at lower energies as the cascades 
redistribute radiation from the higher energies to the lower 
energies (see e.g., Strong et al. 1973, 1974; Coppi & Aha- 
ronian 1997; Kneiske & Mannheim 2008; Inoue & Totani 
2009; Venters 2010). In considering the spectral subpopula- 
tions of blazars, flat-spectrum radio quasars (FSRQs) and 
BL Lacertae-type objects (BL Lacs), the intrinsic spectra 
of FSRQs are much softer than that of both BL Lacs and 
EM cascades, and their relative prominence with respect to 
the EGB can change with energy (Venters 2010; Venters 
& Pavlidou 2011). As such, the fluctuation angular power 
in the EGB should change with energy even if FSRQs and 
BL Lacs have similar anisotropy properties. Furthermore, as 
noted earlier, the deflection of EM cascades by the IGMF 
can spatially extend the emission around gamma-ray sources 
and reduce the observable cascade radiation. Thus, the im- 
pact of a strong IGMF is to reduce the blazar anisotropy 
both because of the increased angular scale of the sources 
and because of the reduction in the overall blazar contri- 
bution to the EGB (assuming that blazars are the most 
anisotropic EGB component). On the other hand, in the 
absence of a strong IGMF, the EM cascades are highly col- 
limated in the direction of the primary photon (angular 
spread ~ m e /E e ), and their impact on the blazar anisotropy 
is to reinforce it due to the enhancement in the overall blazar 
contribution to the EGB. Hence, in the two extreme cases 
for the IGMF, EM cascades could have a profound and oppo- 
site impact on the blazar anisotropy and, by extension, the 
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EGB anisotropy. Therefore, studies of the EGB anisotropy 
could prove a sensitive probe into the nature of the IGMF. 

In this paper, we explore the impact of EM cascade ra- 
diation from blazars on the EGB anisotropy in the extreme 
cases of a negligible IGMF and an IGMF that is sufficiently 
strong to essentially isotropize the cascade emission (i.e., 
from the observationally practical standpoint). This study 
is meant to serve as a proof of concept for the feasibility of 
using EGB anisotropies as an IGMF probe. In Section 2, we 
present the methodology and parameters employed to deter- 
mine the anisotropy energy spectrum and the EM cascades. 
In Section 3, we present the results of the calculation. 


2 METHODOLOGY 

2.1 Collective Intensity of Unresolved 
Astrophysical Sources 


The contribution to the EGB due to unresolved blazars can 
be viewed as the superposition of the collective intensity of 
intrinsic blazar spectra and the contribution from the cas- 
cade radiation from the interactions of VHE photons from 
blazars with the EBL: 


Ie(Eo) = r E (E 0 ) + I c e (E 0 ) , 


(1) 


where the intensity, Ie(Eo), is given in units of photons per 
unit area per unit time per unit solid angle per unit energy 
emitted at observer frame energy Eq. In calculating the col- 
lective spectrum of unresolved blazars, we follow the formal- 
ism outlined in Venters & Pavlidou (2011). Blazar spectra 
are taken to be smoothly broken power laws: 


Fe(Eq) = Fq 



(2) 


where Fe(Eq) is the differential photon flux in units of pho- 
tons per unit area per unit energy per unit time, E b is 
the break energy, ai = a — Aai is the low-energy slope, 
ai = Q + Aa 2 is the high-energy slope, a is the gamma-ray 
photon spectral index of the blazar assuming an unbroken 
power-law spectrum, 2 and n quantifies the sharpness of the 
transition from the low-energy power law to the high-energy 
power law (taken to be 1). 

The total flux, F, of photons with energies greater than 
some fiducial energy, Ef (taken to be 100 MeV), is found by 
integrating Fe(Eq) over energy, 


F = F 0 


f 


EoV 1 (Eo 
E b \E b 


dEo . 


( 3 ) 


Then, the contribution of a single unresolved blazar to the 
EGB is 

7 = F[(E 0 /E b ) ai + (Ep/Eb) 012 ] -1 

1 4t r/“ [(Eo/Ebpi + {Eo/Ebpi}- 1 dE 0 


where the flux of one source is uniformly distributed over 
the entire sky in anticipation of an isotropically distributed 


cosmological population. Including absorption due to the 
EBL, the contribution becomes 

7 = F [(Eo/E b )^ + (Ep/Et)^]- 1 r lEo ,*) (5] 

1 47 r/“ [(Eo/Ebpi + {Eo/Ebpi}- 1 dE 0 

where r(Eo,z) is the optical depth for a photon with ob- 
server frame energy Eo, emitted at redshift 2 . In determin- 
ing the total unresolved blazar contribution including EBL 
absorption, one would generally rewrite Equation 5 in terms 
of the blazar gamma-ray luminosity and redshift, multiply 
by the blazar GLF and intrinsic spectral index distribution 
(corrected for the spectral bias endemic to flux-limited sur- 
veys 3 as per Venters et al. (2009)), and finally integrate over 
gamma-ray luminosity, redshift, and spectral index, thereby 
accounting for evolution in both the sources and the EBL 
(through the redshift-dependent optical depth; see e.g., Ven- 
ters et al. 2009). However, in performing a calculation that 
includes cascades, the population of blazars must be divided 
into its subpopulations of FSRQs and BL Lacs since these 
subpopulations have distinct spectral properties. In order 
to perform a full calculation that accounts for the evolution 
of sources (as in Venters (2010)), we would require reliable 
GLFs for each subpopulation. In the case of FSRQs, a reli- 
able GLF does exist (Ajello et al. 2012b), but one has yet 
to be determined for BL Lacs since it is difficult to obtain 
reliable redshifts for them. While BL Lac GLFs based on 
EGRET and early Fermi data do exist (see e.g., Inoue & 
Totani 2009), a BL Lac GLF based on the current and more 
sensitive Fermi dataset remains in the works (Ajello et al. 
2012a), and it is as yet unclear how the possible modifica- 
tion of BL Lac spectra by cascades could impact the inter- 
pretation of a LAT-measured BL Lac GLF (though, work 
on this topic is currently being performed by the authors). 
Given that the purpose of this paper is to serve as a proof- 
of-concept for using relative changes in the EGB anisotropy 
to constrain the IGMF, for the sake of simplicity, we do not 
incorporate full blazar GLFs in our calculations. Rather, we 
assume typical luminosities and redshifts for all blazars, de- 
ferring a more detailed calculation to the future, at which 
time reliable GLFs for both FSRQs and BL Lacs will have 
been measured. As such, the only integral that must be per- 
formed is that over the spectral index: 


Ie(Eo) = 



[{Ep/Eb)^ FjEo/Eb )^]- 1 

S(E f ,a) 


(6) 


where Io is a normalization constant, p(a) is the measured 
spectral index distribution (SID), and 



We fix the normalization constant so that the total astro- 
physical (blazars, SF galaxies, and EM cascades) contribu- 
tion to the EGB fits the spectrum of the EGB as measured 
by the Fermi - LAT (Abdo et al. 2010d). 


2 For the sake of clarity, we use a for the photon index rather 
than r, which is commonly used in the literature. 


3 We note that the Fermi survey is not exactly a flux-limited 
survey due to the non-uniformity of the total diffuse background 
throughout the sky. However, for the purposes of this paper, we 
neglect this effect. 
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In determining the measured SID, we make use of the 
results presented in Venters & Pavlidou (2011) from a likeli- 
hood analysis fitting Fermi-LAT First Catalog (Abdo et al. 
2010f) FSRQs and BL Lacs to Gaussian SIDs accounting for 
errors in measurement of individual blazar spectral indices. 
For these blazar subpopulations, the maximum-likelihood 
Gaussian SIDs can be characterized by means (ao) and 
spreads ( 00 ) with parameters determined to be ao = 2.45 
and ao = 0.16 for FSRQs and ao = 2.17 and ao = 0.23 for 
BL Lacs. Following Venters & Pavlidou (2011), we model 
the spectral breaks by taking Aon = 0.1, Aa 2 = 0.9, and 
Eb,o = 4 GeV for FSRQs, though we take Eb,o = 10 TeV for 
BL Lacs (treating them as High-frequency-peaked BL Lacs, 
which do not exhibit breaks in the LAT energy range; see 
e.g., Abdo et al. 2010c). 

In principle, one would prefer to perform the above pro- 
cedure on subsamples of local universe blazars to ensure that 
only intrinsic spectra are used (i.e. , that the spectra have 
not already been modified by EBL absorption and EM cas- 
cades). However, the population of local universe sources is 
sparse, and the results of a likelihood analysis performed 
on such a population would be subject to a considerable 
degree of uncertainty. Even so, in the case of FSRQs, the 
spectra are too soft to produce an appreciable amount cas- 
cade radiation (Venters 2010). For this reason and because 
the photon index of the spectrum of cascades is relatively 
hard (F ~ 2; see Venters 2010 and Section 3 of this paper), 
we do not expect substantial modifications of FSRQ spectra 
due to EM cascades (conversely, we do not expect a blazar 
with a hard measured spectrum to have a much softer intrin- 
sic spectrum). Notably, the spectra of local universe FSRQs 
appear to be similar to that of the FSRQ population at 
large (Abdo et al. 2010f). On the other hand, the spectra of 
BL Lacs are sufficiently hard that one might be concerned 
that the spectra are substantially modified by EM cascades. 
Nevertheless, we note that the population of local universe 
BL Lacs appear, on average, to have harder spectra than 
the BL Lac population at large. As such, rather than at- 
tempting a likelihood analysis on the sparse local universe 
BL Lac population (which also have large measurement un- 
certainties in their spectral indices), we make use of the SID 
for the whole BL Lac population as this represents a more 
conservative choice. 

Finally, while the use of broken power-law spectra is 
motivated by our fits to the gamma-ray spectra of other 
types of blazars in the LAT energy range (Venters & Pavli- 
dou 2011), we should note that typical one-zone Synchrotron 
Self-Compton models of BL Lacs predict their gamma-ray 
spectra to be more curved (e.g., log-parabola). If the gamma- 
ray spectra of BL Lacs are more curved than the spectra se- 
lected here, one might expect there to be less cascade than 
predicted here since there would be fewer high-energy pho- 
tons. On the other hand, the high-energy index beyond the 


break is already pretty soft (~ 3) in our template spectra, so 
we do not expect the amount of cascade to change substan- 
tially if we selected log-parabola rather than broken power- 
laws. However, we acknowledge that the overall choice of 
template spectra (e.g., low-energy and high-energy indices 
and break energy) could affect the amount of cascade emis- 
sion, and its impact on our results is currently being studied 
in more detail. 

Since the AGN population peaks ~ 1, we take the typ- 
ical redshift of blazars to be z B rc = 1. In modeling the sup- 
pression of the source brightness due to magnetic deflection 
of cascades (see Appendix), we take the jet opening angle of 
blazars to be ~ 1°. The model of the EBL was taken from 
Franceschini et al. (2008). 

We also include a model of the SF galaxy contribution 
to the EGB. In this case, we assume that the gamma-ray 
luminosity of an SF galaxy is mainly due to the decay of 
neutral pions created through galactic cosmic-ray interac- 
tions with neutral hydrogen in the galaxy. It is thought that 
galactic cosmic rays are accelerated in supernova remnants 
and then diffuse throughout the galaxy. Hence, the flux of 
galactic cosmic rays is expected to be proportional to the 
galaxy’s supernova rate, which is, in turn, proportional to 
its star-formation rate. Furthermore, the neutral hydrogen 
gas provides the fuel for forming stars. As such, the expec- 
tation is that an SF galaxy’s gamma-ray luminosity can be 
parameterized in terms of its star-formation rate (Pavlidou 
& Fields 2002; Abdo et al. 2010a; Fields et al. 2010; Makiya 
et al. 2011; Stecker & Venters 2011; Ackermann et al. 2012b). 
For the purposes of this paper, we model the SF galaxy con- 
tribution using the infrared luminosity function model of 
Stecker & Venters (2011). In this model, the star-formation 
rate of an SF galaxy is assumed to be proportional to its in- 
frared luminosity; thus, the gamma-ray luminosity of an SF 
galaxy is proportional to a power of its infrared luminosity. 
In Stecker & Venters (2011), the relationship between the 
gamma-ray luminosity of an SF galaxy and its infrared lu- 
minosity was determined by fitting the measured gamma-ray 
luminosities of Fermi-detected SF galaxies to their infrared 
luminosities. The SF galaxy contribution could then be de- 
termined by integrating over an infrared luminosity function 
of SF galaxies, which they took from Hopkins et al. (2010). 

2.2 Electromagnetic Cascades from Blazars 

The cascade intensity is given by 
d A N c 

T c _ U iV fo\ 

E dtdAdQdE 5 v ’ 

where dN c /dE is the spectrum of cascade photons due to 
pair production and IC scattering (for full equations and 
derivations, see Venters 2010): 


dN c 

dEo 


Fornax /’^'p,max P 1 .j2 iu 


dN ri (E 0 (l + z)) dNr 2 (E 0 ( 1 + z)) 


dE 


dE 


e -r{E 0 ,*) dfdEpdZj 


( 9 ) 
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where dNr /dE is the spectrum of IC scattered radiation per 
electron of Lorentz factor, T, P(f\E p ,z) is the probability 
that the pair production interaction at a given redshift, z, 
of a primary photon of energy E p will produce electron-type 
particles of energies E e i = / x E p and E e 2 = (1 — /) x E p , 
Fi = / x E p /mc 2 , r 2 = (1 — /) x E p /mc 2 , and d 2 V 7 /dzdE p 
is the continuous spectrum of photons undergoing pair pro- 
duction interactions as a function of redshift and primary 
energy. 

The full numerical integration of Equation 9 is compu- 
tationally cumbersome. As such, in order to calculate the 


cascade spectrum, one typically makes use of a Monte Carlo 
code, such as Cascata (Venters 2010), that propagates pho- 
tons from known spectra and redshift distributions and cal- 
culates the cascade spectrum of each photon. However, for 
the purposes of this paper, we employ a more semi-analytical 
model of the cascade development. Since we have assumed a 
typical redshift for all blazars, we have eliminated the need 
for randomly generated source redshifts. Instead, the spec- 
trum of pair-producing photons, d 2 N^/dzdE p , can readily 
be determined from the absorbed collective intensity of un- 
resolved blazars binned in energy: 


/■Bo, max , rBo.max f 00 !/ CV / p \ai 

I lf s (Eo)dE 0 = I 0 I I dap(a) [{Eo/Eb) 

J En min J En m j r , J — OO 


+ (E 0 /E b ) c 
S{E f ,a) 


1 - e 


-t(E 0 ,z) 


dEo . 


(10) 


For each energy bin, we determined a characteristic red- 
shift for the start of the cascades based on the criterion 
that the characteristic redshift is that at which the num- 
ber of photons that have been “absorbed” in propagating 
from the source redshift is half the number of the total pho- 
tons “absorbed” in propagating to the present epoch. Using 
this criterion, the characteristic redshift can be calculated 
numerically from 


r*(E 0 , z *) = In 


1 

2 


^g r (Bo jZsrc) 



( 11 ) 


Also, we assume that each secondary electron produced in a 
pair-production interaction carries half of the energy of the 
primary photon. Finally, since the number density of CMB 
photons is much higher than those of the other wavelengths 
that comprise the EBL, we assume that the seed photons for 
the IC interactions are CMB photons. In so doing, we can 
neglect Klein-Nishina effects in determining the spectrum of 
scattered radiation. The spectrum of scattered radiation per 
electron of Lorentz factor, T is given by (Venters 2010) 

dNy _ r d 3 N r i _ f d 3 N r , 
de\ J dtdeide J dtdeide 

where T' is the electron Lorentz factor at time, t, t is the 
initial soft photon energy, ei is the scattered photon energy, 
d 3 N r / / dtdeide is the differential scattered photon spectrum 
given by (Blumenthal & Gould 1970) 


dt 


dE e 


dE e de , 


(12) 


d 3 N r , 
dtdeide 


3otc n(e) 

16F' 4 A? 


2ei In 


2 

| | A W2 ^1 

4r ,2 e + £l + £ ~ 2T© 


(13) 


n(e)de is the number density of soft photons, and \dE e /dt\ 
is the rate of energy loss of the electron given by 

dE e 4 /2 

~df = ~ 3 aTcr / m ( e ) de ’ (14) 

We note that this formalism assumes continuous energy 
losses for the electrons, which is appropriate for the regime 
that we are considering. Energy losses due to synchrotron 
emission are also included, but are negligible compared with 
IC losses for the magnetic field strengths considered. Fur- 
thermore, this simplified approach has the advantage that 
recycling cascades (in which IC photons up-scattered to high 


energies are also allowed to initiate EM cascades) is more 
feasible than in the full Monte Carlo approach. As such, for 
BL Lacs, we also calculated the spectrum of second gener- 
ation cascades. Notably, the effect of EM cascades is to re- 
distribute high-energy radiation to lower energies; as such, 
the high-energy component of second generation cascades 
necessary to produce higher generation cascades is greatly 
reduced. Thus, the contribution from generations beyond 
the second generation is negligible. For FSRQs, the intrinsic 
spectra of the sources are sufficiently soft that even the first 
generation of cascades makes a small contribution to their 
overall collective intensity. 4 The suppression of the source 
brightness due to magnetic deflection of cascades (see Ap- 
pendix) is also calculated, and we assume a uniform mag- 
netic field that evolves solely due to the expansion of the 
Universe [B(z) ~ -Bo(l + z) 2 - Neronov & Semikoz 2009], 
where Bo ~ 5 x 10 -14 G. We should also note that for the 
purposes of this paper, we assume that the correlation length 
of the magnetic field is much larger than electron cooling 
distances. Alternative scenarios for the magnetic field and 
source spectra will be explored in future publications. 


2.2.1 A Note About Plasma Instabilities 

Regarding the stability of EM cascades during propagation, 
Broderick et al. (2012) present the very interesting notion 
that plasma interactions between cascades and the inter- 
galactic medium would introduce energy losses that are more 
significant than IC energy losses, thereby quenching the cas- 
cades. Unfortunately, this idea was premised on the asser- 
tion that cascades constitute plasma beams, even though 
the connection with plasmas is not obvious despite cascades 
containing charged particles. In order for a collection of par- 
ticles to constitute a plasma, the particles must be able 
to influence each other electromagnetically (for this reason, 
the particles in a plasma must necessarily be charged parti- 
cles). From our understanding of relativistic beams of parti- 


4 We note that the small amount of cascade radiation is not sim- 
ply the result of the low break energy chosen for this analysis; 
Venters (2010) modeled FSRQ spectra as unbroken power laws 
up to ~ 100 TeV and also found that the FSRQ cascade contri- 
bution was small. 
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cles, the electromagnetic interactions between particles are 
largely suppressed at very high energies (Chao & Tigner 
1999). We can see this by considering the held around a sin- 
gle charged particle at rest and then boosting in the direc- 
tion of propagation. For large Lorentz factors (such as those 
involved in the EM cascade calculations), length contrac- 
tion compresses the electric and associated magnetic Helds 
so that they exist only in a plane perpendicular to the direc- 
tion of propagation of the particle. Since the bulk of the par- 
ticles in the cascades will not lie in this plane, they will not 
be electromagnetically connected to the particle in question. 
For the few particles that could lie in this plane, the Lorentz 
forces exerted on them largely cancel because they are rel- 
ativistic and propagate co-linearly resulting in near cance- 
lation between the electric and magnetic Held components. 
For simple charge distributions (e.g., the two-dimensional 
Gaussian distribution), analytic expressions have been de- 
rived for the Lorentz forces on particles in a beam due to 
other particles within the beam, and such forces have been 
demonstrated to be suppressed by a factor of y -2 (Chao & 
Tigner 1999). As such, even for electrons with energies ~ 
few MeV, the strengths of the forces exerted on them are 
suppressed to the percent level; at energies ~ hundreds of 
MeV and above, these forces are negligible 5 . 

Another method for assessing the impact of collective 
behavior on pair cascades begins with calculating the rela- 
tivistic Debye length for a pair cascade, 


f -y 2 k B T cas \ 1/2 
\ 4-7rne 2 J 


(15) 


where T cas is the transverse temperature of the cascade as 
measured in the laboratory frame (the frame of the inter- 
galactic medium) and n is the particle number density of 
the cascade. The expression for the cascade transverse tem- 
perature is given by (Broderick et al. 2012) 


t' \ J , Teas 
m e c 2 


p\\ fp±\ 2 

2 m e c 2 \pn J 


~5x 10~ 7 



(16) 


where py and p± are the longitudinal and transverse mo- 
menta, respectively, and E 7 is the energy of the primary 
photon. For the particle number density, we use the expres- 
sion given by Miniati & Elyiv (2012), 


n ~ 3 x 10 25 cm 3 


{ E-jL-y \ 

\ 10 45 ergs/s ) 



(i +*r 4 , 


(17) 


for which we take the mean free path for pair production 
for a TeV photon to be 0.8 Gpc at 2 = 1 rather than using 
the full expression. Taking the reference values and plugging 
the results of Eqns. 16 and 17 into Eqn. 15, we find /Debye 
to be ~ 1 kpc. For comparison, we can calculate the radius 
of the cascade, a ~ 8D e , where 8 ~ m e /E e ~ 7 is the 
angular dispersion of the cascade in the absence of magnetic 
fields, D e is the electron cooling distance assuming that IC 
scattering of CMB photons is the dominant cooling process 


(Neronov & Semikoz 2009) 6 , 


3 m 2 c 3 

4(7 tFcmb E e 


10 23 (l + 2) -4 


E e 

10 TeV 


-1 

cm, (18) 


and we assume the electron energy is about half of the pri- 
mary photon energy, E e = E^/2. Assuming a primary pho- 
ton energy of 1 TeV and 2 = 1, we find the radius of the 
cascade to be ~ 0.04 pc. Thus, the Debye length of the 
cascade is much larger than its radius, so charge screening 
is ineffective, and the dynamics of the cascade are domi- 
nated by single-particle behavior rather than collective ef- 
fects (Reiser 1994). Since the mean interaction rate for IC 
scattering of CMB photons is much larger than that for, say, 
bremsstrahlung interactions with protons in the intergalac- 
tic medium (Protheroe et al. 1995), we conclude that the 
process that dominates the cooling of pairs in cascades is 
indeed IC scattering. 

In the case of the nonzero IGMF, the compression of the 
electric and magnetic fields to a plane still applies. Due to 
magnetic deflection, the density of charges falls, suppressing 
the likelihood of finding nearby charges and further increas- 
ing the Debye length. Therefore, we conclude that cascades 
do not constitute the necessary collective phenomena for 
substantial plasma interactions. 


2.3 Anisotropy Energy Spectrum 

The anisotropy energy spectrum, Ci(E), is defined as the 
intensity fluctuation angular power (in units of sr) at a given 
angular scale, t, as a function of energy. In the simple case 
of a two-component background, the total angular power is 
given by (Siegal-Gaskins & Pavlidou 2009) 

C\ ot = ftcY + /fcf 5 + 2 / 1 / 2 C'j lx2) , (19) 

where C Y is the angular power spectrum of component 
(n), f n = I n (E)/Itot(E) is the energy-dependent fraction of 
the total emission arising from component (n), and (j( nxrn l 
is the cross-correlation term for components (n) and (m). 
As CY is a measure of angular fluctuations in units of the 
mean, it is independent of energy for a single population 
of sources with identical spectra. If the two components are 
uncorrelated, C^ nxm ^ = 0, and Equation (19) is further sim- 
plified, 

Cl ot = flC { Y + flC (2) . (20) 

If additional components of negligible angular power con- 
tribute to the total background signal, then these will not 
result in additional terms in Equation (20), but rather they 
will affect Cj ot (E) by changing the fractional contributions, 
/, of each of the components that do contribute to the an- 
gular power. 

In the simple scenarios we explore, we always consider 
two components that do contribute to the total angular 
power: BL Lacs and FSRQs. Because in our adopted model 
the two contributions do not add up to the total EGB sig- 
nal measured by Fermi (Abdo et al. 2010d), we assume 


5 For this reason, it has been noted that collective effects are also 
negligible for high-energy charged particle beams such as those 
made in the Large Hadron Collider at CERN. 


6 Note that since the bulk of the cascade is produced in the first 
generation, we consider only one generation for this simple calcu- 
lation. 
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that the rest of the EGB signal originates from a very low- 
Ct population, such as SF galaxies. For simplicity, we take 
Cf F > = 0 for this component; hence, their contribution only 
enters Equation (20) through the total background intensity, 
Itot(E) (i.e. , the denominators of /i and f%)- As discussed 
in the previous sections, the energy spectrum of the pri- 
mary emission from FSRQs is very soft at very high energies, 
so the cascade signal from FSRQs is negligible; all cascade 
emission that we consider comes from BL Lac VHE photons. 

In the limit of low B , there is little angular spread in the 
cascade photons, so these photons appear to originate from 
the same sources as the primary emission. In this case, the 
only effect of the cascades is to alter the energy spectra of 
blazars, leaving their anisotropy properties unchanged. As 
such, the fractional components of Equation (20) become 
fbi(E) = [I^(E) + I^ C (E)] /Itot(E), while Cj bl) remains 
unchanged. 

In the limit of strong B (Bo ~ few x 10~ 14 G), we 
determined the scattering angles of the cascades (see Ap- 
pendix) and found that for observed photon energies less 
than ~ 100 GeV, the scattering angles were sufficiently large 
such that within the angular spread of the cascade emission 
on the sky (halo radius ~ several degrees at E 1 ~ tens of 
GeV), there would be enough sources to result in overlapping 
(confused) halos 7 . Under these conditions, the cascade emis- 
sion from the entire population of BL Lacs would be fairly 
isotropic, so we take the cascade signal as having angular 
power, Cf asc) = 0. Since the cascades now have an angular 
power that is different from the blazars and the null IGMF 
case, they will enter Equation (20) in the same manner as 
the SF galaxies, only through the total background inten- 
sity, Itot(E). The numerators of /i and fa in Equation (20) 
will then only include the intrinsic emission from FSRQs 
and BL Lacs, respectively (i.e., fb\(E) = 1^7 (E) / Itot(E)) . 
We emphasize that it is for this reason that the total angular 
anisotropy energy spectra of the two magnetic held scenarios 
will be different. Even if the overall EGB cascade contribu- 
tions in the two scenarios were the same, the total angular 
anisotropy energy spectra would be distinct, depending on 
the prominence of the cascades in relation to the other con- 
tributions to the EGB. 

Since, as stated in Section 2.1, a reliable BL Lac GLF 
has yet to be constructed, a full calculation of the angu- 
lar power of BL Lacs as per Ando et al. (2007) is not 
yet feasible. Instead, we normalize our anisotropy energy 
spectra such that in both scenarios, C( 0t ~ 10~ 5 sr around 
1 GeV, consistent with the Fermi measurement of angu- 
lar anisotropies in the EGB at these energies (Ackermann 
et al. 2012a). This results in a blazar angular power 8 of 
£,(fsrq) _ ^BLLac) _ g x ^g-s srj which w iU be testable in 
future publications once a reliable BL Lac GLF is found. The 
measurement at 1 GeV was chosen to balance the competing 
needs of adequate angular resolution and sufficient photon 
statistics, and to minimize possible contamination from any 
remaining galactic foreground emission. The uncertainty at 


7 As based on the number of sources above the Fermi - LAT sen- 
sitivity above 10 GeV. Since the blazar contribution to the EGB 

is the collective emission of unresolved blazars, this criterion is 
quite conservative. 


1 GeV is relatively small, though more data from Fermi will 
likely improve it. 


3 RESULTS 

In Figure 1(a), we plot the spectra of the modeled contribu- 
tions to the EGB. The “total” contribution is the sum of the 
contributions from blazars, SF galaxies, and one of the cas- 
cade models (the total is similar for both models since the 
largest differences in the cascade contributions in the two 
scenarios occur at energies at which the EGB is dominated 
by blazars and SF galaxies). The blazar component hard- 
ens above ~ few GeV due to the increased contribution of 
BL Lacs with respect to FSRQs, as BL Lac spectra remain 
hard while FSRQs soften. The cascade spectrum in the case 
of the zero IGMF is, as expected, quite flat (dN/dE oc E~ 2 ; 
see e.g., Strong et al. 1974; Venters 2010) demonstrating the 
cascade effect of redistributing radiation from high energies 
to lower energies. The sum of the three components fits quite 
well the spectrum of the EGB as measured by Fermi. No- 
tably, the spectra of the cascades for the two cases of the 
IGMF, while similar at the highest energies, are quite dis- 
similar at the lower energies. This is due to the loss of some 
cascade radiation in the strong IGMF case as significant 
amounts of the cascades are deflected away from the line of 
sight of the observer. However, it should be noted that given 
our simple model of the blazar contribution to the EGB, we 
were only able to account for cascade losses due to deflec- 
tion. That is, in our simple picture, we are calculating the 
amount of radiation that an observer looking head-on at a 
given source would be able to detect. It is likely that just 
as cascades from one aligned source are deflected out of the 
observer’s line of sight, cascades from other aligned sources 
and even misaligned sources are deflected into the observer’s 
line of sight. Nevertheless, we note that even if the strong 
IGMF cascade radiation were comparable to that of the zero 
IGMF case, the angular power in the two cases would remain 
distinct since the angular properties of the cascades are dif- 
ferent (see Section 2.3). 

In Figure 1(b), we plot the anisotropy energy spectra 
of the EGB for the two cases of the IGMF. Notably, there 
is an appreciable difference between the two cases starting 
at ~ few GeV. This is due to the fact that above these en- 
ergies, the contribution from EM cascades relative to the 
other components is much more significant than at lower 
energies. As a result, the impact of the isotropization of the 
cascades by the IGMF is much more significant at these 
energies. Aiding in the distinction is that the cascades aug- 
ment the anisotropy in the case of zero IGMF, because they 
enhance the contribution of the most anisotropic compo- 
nent. For comparison, the Fermi - LAT measurement of the 
anisotropy energy spectrum of the EGB (Ackermann et al. 
2012a) is also plotted. With just two years of data, the Fermi 
measurements seem to already favor the higher IGMF case 
over the null IGMF case. 


8 Note that for a Poisson distribution of sources, the angular 
power at a given energy is constant with multipole. 
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Figure 1. a: The spectrum of the modeled EGB (thick solid line) together with individual components: intrinsic blazar emission, 
including both FSRQs and BL Lacs (thin solid line); SF galaxies (dashed line); electromagnetic cascades for zero IGMF (dot-dashed line) 
and nonzero IGMF (double dot-dashed). For reference, the spectra of the EGB based on both Fermi (filled circles; Abdo et al. 2010d) 
and EGRET data (open squares; Sreekumar et al. 1998; open triangles; Strong et al. 2004) are also plotted, b: The anisotropy energy 
spectra of the modeled contributions to the EGB. The solid line is the total model assuming zero IGMF. The double dot-dashed line is 
the total model assuming nonzero IGMF. Grey boxes represent the Fermi - LAT measurement of the anisotropy energy spectrum of the 
EGB (Ackermann et al. 2012a). 


4 CONCLUSIONS & DISCUSSION 

We have studied the effect of a strong IGMF on EM cas- 
cades from blazars and the resulting impact on the cascade 
component of the EGB and the anisotropy energy spectrum 
of the EGB. We have shown that while the spectrum of 
cascade radiation in the case of zero IGMF is flat in E 2 , 
the spectrum in the case of a strong IGMF is suppressed at 
low energies, though there could be an additional contribu- 
tion arising from cascades from other sources, such as radio 
galaxies and propagating UHECRs. 

For the two extreme cases of IGMF strength we have 
considered here, we have calculated the anisotropy energy 
spectrum of the EGB. In the case of zero IGMF, the ef- 
fect of cascades is to augment the anisotropy arising from 
blazars since the cascades increase the contribution to the 
EGB arising from blazars and are highly collimated with 
the primary emission. In the case of the strong IGMF, the 
cascades are deflected substantially and the angular proper- 
ties of the parent blazar population are blurred. Thus, in the 
case of a strong IGMF, the effect of cascades is to reduce the 
anisotropy of the EGB as they hinder the anisotropy arising 
from blazars rather than reinforcing it. At energies at which 
the relative contribution to the EGB arising from EM cas- 
cades is significant, the anisotropy energy spectra of the two 
cases are distinguishable, and current Fermi data already 
seem to favor non-negligible IGMF values. With more data 
and further refinement of this technique, Fermi will be able 


to provide tighter constraints and possibly even measure the 
IGMF. 

It is worth noting that if blazars and their correspond- 
ing cascades make a large contribution to the Fermi EGB 
spectrum above 1 GeV, then the null IGMF scenario would 
be in tension with current Fermi anisotropy measurements 
(see Figure 1(b)). Similar observations have been made by 
Cuoco et al. (2012) and Harding & Abazajian (2012), lead- 
ing both groups to conclude that blazars can explain no more 
than 20% of the EGB above 1 GeV lest they overproduce the 
EGB anisotropy. However, our results demonstrate that in 
the scenario of a non-zero IGMF, it is possible that blazars 
and their corresponding cascades can contribute substan- 
tially to the Fermi EGB spectrum above 1 GeV and even 
reproduce the EGB spectrum above 10 GeV without over- 
producing the EGB anisotropy. 

Many possible contributions to the EGB have been 
identified in the literature, and a multi-component back- 
ground is plausible in light of current Fermi data. The ex- 
tent to which other contributions will impact the differ- 
ence in anisotropy energy spectra in the two cases consid- 
ered here depends on their angular power and whether they 
also have appreciable cascade contributions. For instance, 
other gamma-ray-emitting radio galaxies for which the jet 
is not aligned with the observer’s line of sight (“misaligned” 
blazars) could also contribute to the EGB (Inoue 2011). 
For those radio galaxies that are too misaligned for the ob- 
server to receive substantial intrinsic gamma-ray emission, 
a nonzero IGMF could deflect cascades into the observer’s 
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field-of-view, giving rise to an additional contribution to the 
EGB, provided that the spectra of these sources are hard 
enough to initiate cascades. These sources will be substan- 
tially fainter than blazars, but are not likely to be as nu- 
merous as SF galaxies. Hence, they are not likely to have 
substantial angular power and will enter the calculation of 
the total angular power of the EGB in a manner similar to 
that of the SF galaxies, solely by virtue of their relative con- 
tribution to the intensity of the EGB. We will return to the 
impact of radio galaxies in a future publication. 

In addition to that from radio galaxies, there could also 
be a contribution to the EGB from the annihilation of dark 
matter in the Galaxy. Even if gamma rays from dark matter 
annihilation do not dominate the EGB in any energy band, 
they could have a substantial impact on the total angular 
power of the EGB provided that they make a significant 
contribution to the EGB intensity. The impact of dark mat- 
ter annihilation on both the intensity and angular power 
of the EGB depends on the mass of the dark matter par- 
ticle and the combined cross section for annihilation into 
channels that result in gamma rays. Some models of the 
EGB contribution from dark matter can result in spectra 
that peak at the same energies at which the EM cascades 
contribute most significantly to the EGB. However, the sig- 
nificance of the dark matter contribution to the EGB also 
depends on the amount of substructure present in the dark 
matter halo with the largest signals occurring in scenarios 
in which the halo contains numerous large clumps (Siegal- 
Gaskins 2008). In such scenarios, the angular power of dark 
matter is substantial and will likely enhance the total angu- 
lar power of the EGB at the energies at which the dark mat- 
ter signal peaks (Siegal-Gaskins & Pavlidou 2009). In both 
of the IGMF cases considered in this paper, the impact of a 
substantial contribution from clumpy dark matter at ener- 
gies of tens of GeV would be to increase the angular power 
in the last energy bin in Figure 1(b). The angular power in 
the two scenarios would still be distinct, though the distinc- 
tion would depend on the significance of the dark matter 
contribution relative to that of blazars and EM cascades. 
Nonetheless, in order to push the two models to yield the 
same angular power in the last energy bin, the contribution 
from dark matter would have to be substantial, and even in 
the case of nonzero IGMF, it could prove difficult to explain 
the flatness in total angular power of the EGB with energy 
as measured by the Fermi - LAT (Ackermann et al. 2012a). 
Even worse is the scenario in which there is no substantial 
cascade contribution, but there is a substantial dark matter 
contribution - without a low-anisotropy contribution to mit- 
igate the total angular power, the angular power in the last 
energy bin would rise substantially. We must also consider 
the possibility that the dark matter halo is smooth rather 
than clumpy. In this case, the angular power of dark matter 
would be minimal, but it could be difficult for dark matter to 
contribute significantly to the EGB. In any case, Fermi has 
already begun to constrain models weakly interacting mas- 
sive particles as dark matter candidates, including tantaliz- 
ing hints of a signal at ~ 130 GeV near the Galactic Center 
(see e.g., Weniger 2012). In the future, the much anticipated 
results from the Large Hadron Collider could provide even 
more clues on the nature of dark matter. Thus, the question 
of the dark matter contribution to the EGB could soon be 
answered. 


Finally, there could also be an appreciable cascade con- 
tribution from the interactions of UHECRs with the EBL. 
The impact of cascades from UHECRs on the total angular 
power of the EGB depends on the still-unknown sources of 
UHECRs, but the most likely scenario is that the difference 
in angular power between the two IGMF cases will still be 
present and may, in fact, be augmented. In order for the sep- 
aration in angular power to disappear with the inclusion of 
UHECR cascades, the UHECR sources would have to have 
very little angular power. Such a scenario might have been 
possible in the top-down models of UHECRs, but such mod- 
els are tightly constrained by the upper limits on the UHE 
photon flux determined by the Pierre Auger Observatory, 
and some such models have already been excluded (Pierre 
Auger Collaboration 2009). As such, the current discussion 
of the origins of UHECRs largely plays out in the astro- 
physical arena, and the classes of astrophysical sources that 
are typically thought to be capable of accelerating UHECRs 
tend to be rare and powerful sources (as are blazars). Hence, 
the scenario in which UHECR sources exhibit very little an- 
gular power is unlikely. Furthermore, for electrons with en- 
ergies greater than those considered here, the IC process 
becomes more inefficient as Klein-Nishina effects become 
important. Therefore, the cascade photons from UHECR 
interactions that would contribute to the EGB are likely 
produced by electrons with energies similar to those found 
in the VHE gamma-ray cascades produced here. Further- 
more, UHECRs are themselves charged particles and would 
have already been deflected by the IGMF prior to initiat- 
ing cascades. Thus, the inclusion of UHECR cascades would 
likely also result in a difference in angular power between the 
two IGMF cases, which may even be more prominent than 
demonstrated in this paper. Given the possibility of a large 
contribution to the EGB from UHECR cascades (e.g., Ahlers 
et al. 2010) and the possible implications for the sources of 
UHECRs, we will devote a future publication to the study 
of the impact of UHECR cascades on the angular power of 
the EGB. 

We note that an advantage of our approach is that 
rather than focusing on a few select sources, our method em- 
ploys information from the whole observable sky. As such, 
our approach could prove capable of providing more sensitive 
constraints on the global properties of the IGMF. Further- 
more, our approach is not limited to uniform magnetic fields 
with large coherent lengths; other field configurations can be 
probed with this technique. In practice, in order to use this 
technique to provide an IGMF measurement, intermediate 
cases of IGMF strength need to be considered as well, and 
the uncertainties entering the calculation of the expected 
anisotropy energy spectra in each case carefully evaluated. 
In the case of intermediate IGMF strengths, a numerical 
simulation of the expected angular power is appropriate as 
the angular anisotropy properties of the cascade radiation 
will be altered with respect to the primary BL Lac emis- 
sion, though not completely suppressed as in the case of the 
strong IGMF. In fact, since the cascade radiation is corre- 
lated with the primary blazar emission, the cross-correlation 
term of Equation (19) would have to be calculated. As such, 
for intermediate IGMF strengths, a simple analytical calcu- 
lation of the cascade effect on the total angular power is not 
possible. Concerning uncertainties, these enter our calcula- 
tion primarily through the model assumptions for each EGB 
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component (which can, in turn, self-consistently determine 
their anisotropy properties), and through the EBL model. 
We will expand our investigation of the practical usage of 
the anisotropy energy spectrum of the EGB to constrain the 
properties of the IGMF in future publications. 

In the meantime, this simple, proof-of-concept calcula- 
tion has already demonstrated that the anisotropy energy 
spectrum can be a powerful tool in constraining the IGMF. 
With a few more years of data from Fermi , we can begin to 
probe the IGMF and uncover more clues into the origins of 
UHECRs and large-scale structure formation. 
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Figure 1 . a: Geometry for the angular deflection of two generations of EM cascades, b: Geometry for the angular deflection of one 
generation of EM cascades. 


APPENDIX A: SUPPRESSION OF THE 
SOURCE BRIGHTNESS DUE TO MAGNETIC 
DEFLECTION OF EM CASCADES 


In calculating the suppression of the source brightness due to 
magnetic deflection of EM cascades, we must determine the 
angular spread of the cascades, 9 STC , as viewed by a funda- 
mental observer at the source redshift, z sr c (Figure Al(a)). 
In order to determine d STC for multiple cascades, we begin 
with considering the problem for one generation of cascades 
(Figure Al(b)). In this case, 0\, the angular spread of the 
cascade as viewed by a fundamental observer at zi, can be 
expressed in terms of quantities evaluated at the redshift 
of interaction of the “primary” photon, Z2, and the angular 
diameter distance from z src = z\ to 22, 1X1(21,22): 


0 i 


d 2 

1X4(21,22) ’ 


(Al) 


where (X is the proper lateral distance the electron has trav- 
eled as measured by a fundamental observer at 22 (see e.g., 
Longair 2008 ), and we neglect Klein-Nishina effects, so the 
electron cooling distance, d|, is not cosmological. By geom- 
etry, d2 is given by 

d2 = I2 sin 82 = 2r2 sin 2 £2 , ( A 2 ) 


where I2 is the straight path of d % , r% is the Larmor radius 
at 22 of an electron with energy, Ef, in a magnetic field of 
strength, B(z 2 ) = B 0 (l + z 2 ) 2 , given by (Neronov & Semikoz 
2009 ) 


r-2 = (3 x 10 28 cm) 


Bin) 

10- 18 G 


10 TeV 


(A 3 ) 


The angular diameter distance from 21 to 22 can be found 
from the angular diameter distance from 22 to 21 through 
the Reciprocity Theorem: 

Da(zi, z 2 ) = i*— ' -1X4(22, 21) , (A 5 ) 

1 ~r Z 2 

where Da(z2, zi) for a flat universe is given by 

Da{z 2 , 21) = * - [Dc( 2 l) - Dc(Z2)] , (A 6 ) 

1 Z\ 

and D c (z) is the comoving distance at 2. 

Given 6 1 in terms of known quantities, we can find 9 src 
for two generations of cascades using an equation similar to 
Equation Al. Returning to Figure Al(a), we note that in 
this case, the primary photon is emitted at 2 src and inter- 
acts at 21 emitting electron-type particles that do not prop- 
agate cosmological distances. One electron-type particle is 
deflected away from the straight-line path of the primary 
particle through an angle, < 5 i, and up-scatters CMB pho- 
tons. The up-scattered photons then propagate until they 
interact at 22, creating more electron-type particles. Since 
we neglect generations of cascades higher than the second 
generation, 9\ is given by the previously developed formal- 
ism. The angle through which the electron at 21 is deflected 
away from the straight-line path of the primary photon, < 5 i , 
is also given by Equation A 4 , but for quantities evaluated 
at 21 . 

Given 9 1 and 5 i, the total angular spread of the cascades 
as viewed by a fundamental observer at 21 is 0\ = 9\ + 81. 
The lateral spread of the cascade as measured by a funda- 
mental observer at 22, D2, is given by Equation Al: 

D 2 = #(1X4(21,22) . (A 7 ) 


and 82 = ( 1 / 2 ) 91 is the angle through which the electron is 
deflected away from the straight-line path of the “primary” 
photon, which is equal to half of the magnetic deflection 
angle of the electron given by (Neronov & Semikoz 2009 ) 


9: 


2 = T = (3x 10” 

r 2 


5 )(l+*a)- 


B 0 


10- 


G 


El 


10 TeV 


(A 4 ) 


Finally, we can again use Equation Al to find the angular 
spread of the cascades as viewed by a fundamental observer 
at the source redshift: 


X >2 

Da(z src, 22 ) 


(A 8 ) 


The flux suppression of the cascade due to magnetic deflec- 
tion can be found from the flux of the cascade in the absence 
of magnetic fields by relating the solid angle of the deflected 


© 2013 RAS, MNRAS 000 , 1-12 


12 T. M. Venters and V. Pavlidou 


radiation (determined from # src ), to that of the unaf- 

fected emission (in this case, the solid angle of the jet), Afi: 

^AQ d = ■ ( A9 ) 

We note that this procedure is only valid for non- 
isotropically emitting sources. Conceptually, there should be 
no flux suppression for an isotropic source because the sym- 


metry of the problem demands that whatever emission that 
is deflected out of the observer’s cone of sight is replaced by 
emission being deflected in. 

This paper has been typeset from a TgX/ lATjTjX file prepared 
by the author. 
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